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«) ' Abstract. The high resolution capabilities of the High Energy Stereoscopic System (HESS) intro- 

I duced a new era in Gamma-Ray Astronomy, and opens a new window on pulsar wind nebula (PWN) 

■ research. A rotationally induced jet (associated with PSR B 1509-58) is resolved for the first time in 
_ I 7-rays, allowing us to trace the particle transport directly, without having the complicating effect of 

, spatially varying field distributions on the synchrotron emissivity. For PWN older or more extended 

■ than Crab (i.e. those with lower field strengths), HESS also reveal the properties of electrons con- 
I tributing to the EUV/soft X-ray synchrotron bands, whereas EUV/soft X-rays suffer from severe 
. interstellar absorption effects. Finally, HESS morphological studies of evovled PWN also allow us 

to directly measure the effects of assymetric reverse shock interactions due to SNR forward shock 
I expansion into the inhomogeneous interstellar medium. 
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§ INTRODUCTION 

. Weiler & Panagia [1] first coined the term "filled-center" or "plerionic" nebulae, fol- 

^! lowing the discovery of galactic sources with the following features: (1) filled center or 

O ' blob-like form; (2) a flat radio spectral index between and -0.3; (3) a well organised 

-t— ' ■ magnetic field; and (4) a high integrated linear polarisation at high radio frequencies. In 

! 1980, Weiler & Panagia [2] added Vela X, part of the Vela SNR, as a member of this 

\ class of sources. The offset of the Vela X PWN relative to the Vela pulsar was puzzling, 

until resolved as discussed below. 

X-ray Astronomy was quite instrumental in detecting and characterising these pulsar 
I wind nebulae (PWN) and we were left with Radio and X-ray Astronomy contributing 

mostly to our knowledge of PWN, other than the Crab Nebula. A few useful contribu- 
tions from Infrared Astronomy were also made, filling the spectral gap between radio 
and X-rays. The radiation mechanism for radio through X-rays is generally believed to 
be synchrotron emission as a result of pulsar injected electrons, with X-rays resulting 
from freshly injected electrons, whereas relic electrons contribute mostly to the radio 
part. The ordered magnetic field is usually toroidal (due to the spinning pulsar), with 
magnetization resulting from the deceleration of the post-shock pulsar wind flow [3]. 
The synchrotron brightness is a result of a convolution of the electron density and the 
magnetic field strength, whereas the same electrons also scatter ambient photons into 
the high energy to ultra high energy y-ray range. Both components should be observed 
to understand the dynamics of the pulsar wind. 

The well-known optical "wisps" in the Crab Nebula are believed to mark the pulsar 
wind shock radius, where the pressure from the unshocked (upstream) pulsar wind 
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balances the (downstream) nebular pressure [3]. The high resolution X-ray images of 
Chandra contributed significantly to this field, showing that the presence of a torus 
and perpendicular jet are common to PWN (see [4] and references therein), with an 
underluminous region on the inner part of the torus, marking the PWN shock radius. Ng 
& Romani [4] managed to derive the geometry of these structures, which establishes the 
orientation of the pulsar spin axis (aligned with the X-ray jet) on the plane of the sky. 
Thus, unified pulsar-PWN models should combine measured radio pulse polarisation 
swing data and high resolution X-ray PWN geometry for internal consistency. 

A complicating factor for evolved PWN is the effect of the reverse shock on the PWN 
for SNR expansion into an inhomogeneous ISM: van der Swaluw et al. [5] showed that 
after the passage of the reverse shock into the PWN, the latter will expand subsonically, 
settling to a radius, which is ~ 25% of the SNR forward shock radius. Blondin et al. 
[6] also showed that an inhomogeneous ISM pressure will result in the reverse shock 
returning to the PWN at different times of its evolutionary history, which can shift the 
position of the PWN, as observed for Vela X and G 18.0-0.7 (the latter is associated with 
PSR B 1 823-1 3 [7] .) Since most massive star formation is taking place in inhomogeneous 
molecular clouds, we expect to see several offset PWN, which can be studied with the 
high resolution imaging capabilities of HESS, since the HESS emitting electron lifetimes 
can be longer than the epoch of PWN shift. 

A review of the HESS telescopes and basic overview of results are given by G. 
Hermann in these proceedings. In this paper we concentrate on the PWN results of HESS 
and what we learn from these results, combined with multiwavelength information. At 
the time of the writing of this paper, a number of PWN has been detected by HESS, but 
we only discuss the PWN of SNR G 0.9-1-0.1, PSRB1509-58, and PSRB 1823-13. 



HESS J 1747-281 = PWN OF SNR G 0.9+0.1 

Unresolved TeV emission from SNR G 0.9-1-0.1 was detected after 50 hours of live-time 
observations [8]. The image in Figure la shows that the TeV source is coincident with 
the PWN, and is smaller than the shell size associated with this composite SNR. This 
source is one of the weakest TeV sources detected so far, with an integral flux above 200 
GeV of F(> 200GeV) = (5.7 ±0.7,t,tf ±1.2,,.,) x lO^^^ cm-h'^ and a photon spectral 
index of Thess — (2.40 ± 0. 1 Igtat =t 0.20,^,) . This spectrum is shown in Figure lb [8]. 

de Jager et al. [9] derived an analytical expression for the TeV 7-ray spectrum in 
terms of the X-ray spectrum, assuming that the X-ray spectrum with energy index ax 
(between and 2.5) extends to the EUV/soft X-ray range, since synchrotron emitting 
electrons producing such soft emission also IC scatter CMBR photons into the TeV 
range for a magnetic field strength in the range ~ 10/iG. The IC 7-ray energy spectrum 
(with ^TeV the 7-ray energy in units of TeV) at earth, written in terms of the observed 
monochromatic X-ray flux Siy^y at 1 keV is then (in units of TeV/cm^/s/TeV) 

fl 4 X 10~^)*^* 

52.7K(£Tev) = 6.6 X IQ-^^ ^—^—^exp [2.2a^ - 0. 126a2] SikevCz.vK^T,^, 



with the first-order Klein-Nishina correction function given by C2.7K = [1 — (2.4 -|- 
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FIGURE 1. (a) Left: The 90 cm radio flux map of G 0.9+0.1 with HESS confidence contours: white 
soHd (dashed) circles: observed (simulated point spread function) 40% and 80% peak brightness; green 
dashed circle: 95% upper limit on HESS source size, (b) Right: The HESS photon spectrum of GO.9+0.1. 



3.9a^ + 0.42a^)(o.03£Tev)^''^]- If IC emission on the 25K galactic (or swept-up) dust 
component (with density U25 in units of eV/cm^) is important, we have to add the 
additional IC TeV flux given by the expression: 

S25k{Et,y) = 0.023 (25/2.7) + 3) ( a, + 4) C/25 

with the first-order Klein-Nishina correction now larger [9] C25K = [1 ~ (2-4 + 3.9ax + 
0.42a^)(0.09£'Tev)^'^^]- This correction function is only valid for C < 0.5 and higher 
order corrections are then necessary for larger energies. 

Milky Way type galaxies show that the IR emission from dust dominates the starlight 
component near the galactic center (R. Tuffs & C. Popescu, personal communication, 
2005), contrary to the assumption in [8], where the starlight was assumed to dominate 
the dust component. This observation also makes the interpretation of the TeV spectrum 
simpler as discussed below: 

XMM Newton observations of GO.9+0.1 showed that the energy index Ux = Tx — ^ 
increases with increasing radius (possibly due to synchrotron cooling), giving a value of 
Ux = 1.4 ±0.2 for an outer shell of radius 1.8 to 3.5 pc [10], which is larger compared 
to the 0.6 pc size of the Crab X-ray nebula. It is reasonable to assume that the field 
strength in this outer region would be smaller compared to the compact core, so that 
most electrons have accumulated here over the pulsar lifetime, resulting in the observed 
IC TeV emission. The 95% upper limit on the HESS source radius (containing 68% of 
the events) is also 3.2 pc [8] (assuming a distance of 8.5 kpc to the galactic center.) The 
HESS source therefore includes the total X-ray PWN and we expect most HESS emitting 
electrons to be concentrated in the region where the field is lowest and volume is largest. 
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FIGURE 2. (a) Left: The Chandra image of the PWN of PSRB1509-58, indicated by the extended 
blue jet and compact torus perpendicular to the jet, centered on the pulsar. The red structure at the top 
is the thermal nebula RCW89, undetected by HESS [13]. (b) Right: The HESS image of the PWN of 
PSRB1509-58, showing a similar jet structure as indicated by the projection from Chandra (red and 
white lines) and the ROSAT PSPC contour overlay (white contours) [12]. 



It also appears to be valid to extrapolate the X-ray spectrum into the EUV range, since 
the spectral break between radio and X-rays is only observed around 10^^ to 10^^ Hz 
[11]. The photon index of TeV emission should therefore also be Thess ~ + 1 = 2.4 
(as observed.) Assuming a monochromatic flux of 5ikeV = 0.0014 cm^/s/keV at 1 KeV 
[10], we calculate a magnetic field strength in the large outer volume of the PWN 
as 10/iG, assuming a 25K photon density of 1 eV/cm^ in this region of the galaxy. 
However, a 2.3 times higher density would give a field strength of ~ 14/iG. 

Thus, the agreement between the X-ray and TeV spectral indices ( within uncertainties ) 
exclude the need to invoke a dominant central starlight component, consistent with 
observations of other Milky Way type galaxies. Furthermore, a more detailed IMHD 
approach is needed to put constraints on the injection parameters of the unseen pulsar. 

HESS J1514-591 = PWN OF PSRB1509-58 

This source was easily detected by HESS, with a significance of 25 a after 22 hours of 
observation [12]. Good statisics allowed the measurement of the TeV spectrum with 
high precision: The monochromatic flux at 1 TeV is (5.7 ±0.2stat ± ^-^-sys) x 10^^^ 
cm~^s^^TeV^^ and the power law spectral index is Thess = (2.27 ±0.03ita? ±0.20^^.5) 
for the range 0.28 to 40 TeV. The power law is unbroken [12]. 

From Figure 2 we see that the TeV morphology follows the X-ray morphology in 
terms of the jet-like structure: This pulsar also shows a torus and jet, with the torus 
clearly resolved by Chandra [13]. Contrary to other compact PWN, the jet dominates 
the toroidal component, and the same is seen in TeV y-rays. This is the first time that an 
astrophysical jet has been resolved in the y-ray domain. The TeV emission extends to a 



radius of 0.3° (or 25 pc) from the pulsar, assuming a distance of 5 kpc to tlie pulsar 

Whereas the jet close to the pulsar reveals a photon index of Fen = 1.6 ± 0.1 (with 
"cn"="compact nebula" derived from feature "C" of [13]), the diffuse emission up to 
0.3 degrees reveal a spectral index of Fen = 2.04 ± 0.01 ("en"="extended nebula", [14]). 
This steepening (by Fen — Fen = 0.44 ± 0. 10) is consistent with synchrotron cooling. The 
TeV emission corresponds to this extended cooled region and the additional steepening 
Thess — Ten = 2.27 — 2.04 = 0.23 can be explained by the addition of a 25K dust 
component by galactic and swept-up dust in the SNR of this pulsar [15, 12]. A discussion 
of the cooling break will be given below. 

Aharonian et at. [12] found that galactic dust emission at a level of 2.3 eV/cnr' could 
explain the observed HESS spectrum, assuming a mean magnetic field strength of 17/iG 
in this extended region. Before the HESS era, Gaensler et al. [13] argued that the spectral 
break due to synchrotron losses should be just below 1 keV (for an 8/iG field), whereas 
the field strength of 17/iG derived from TeV observations predicts a spectral break near 3 
TeV, which is not observed [12]. Note however that this cooling break scales as b~^T~^ 
in the Thomson limit, and that the nebular field strength was larger during earlier epochs 
(i.e. expansion causes the field strength to decrease with time). Just a factor two larger 
B (without changing the spindown age), during earlier epochs, will already push the 
spectral break down by a factor of 16. The result is then an unbroken power law down 
to at least 0.3 TeV, as observed. Thus, more realistic PWN models should allow for field 
evolution with time, whereas constant field models may predict a TeV spectral break, 
which is not observed. 

HESS J1825-137 =(?) G 18.0-0.7 = PWN OF PSRB1823-13 

The source HESS Jl 825- 137 was discovered as part of the HESS galactic plane survey 
[16]. It is resolved with an RMS radius of ~ 10', and the 101 ms period Vela- like pulsar, 
PSRB1823-13, is the only candidate source within the area resolved by HESS. The 
HESS Collaboration recently made a case for the association of HESS J 1825- 137 with 
G 18.0-0.7, based on morphological and spectral considerations [17]. 

Despite the fact that the X-ray PWN G 18.0-0.9 is located south of this 21 kyr old 
pulsar, the association between the pulsar and the X-ray PWN is considered to be firm 
[7]. Figure 3a shows an excess slice plot taken along the north-south direction, showing 
this assymetry. Gaensler et al. [7] invoked the Vela X type explanation of Blondin et 
al. [6] to explain the southwards offset of the PWN relative to PSRB1823-13. From 
Figure 3a we can see that the bright X-ray compact nebula (containing freshly injected, 
uncooled electrons) is symmetric with respect to the pulsar, whereas the older (cooled) 
electrons result in a southward offset X-ray nebula. The TeV nebula is also shifted to the 
south (Figure 3b), with an extension (~ 0.5°), which is even larger compared to X-rays. 
This is expected based on the longer synchrotron lifetimes of TeV vs X-ray emitting 
electrons. 

The 0.5° diameter of the TeV source translates to a radius of RvwN ~ 17<i4 pc 
for a distance of J = 4 kpc. The SNR forward shock radius is then expected to be 
RsNR ~ ^RvwN = 70(i4 pc [5], which can be reached if we invoke a relatively large 
initial spindown power and pulsar braking index less than the canonical value of n = 3, 
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FIGURE 3. (a) Left: The XMM north-south excess sHce of the PWN of PSR B 1 823- 1 3 (from [7]), (b) 
Right: similar north-south slice of HESS J 1825- 137, also showing emission south of the pulsar, but over 
a much larger size (from [17].) 



to increase the age of the SNR. The Sedov-Taylor SNR size would also increase with 
increasing age, giving a radius of 



RsNR = VSpc 



EsN 0.003 cm-3\" V i ^ ^-^ 



lO^^Qrgs N J \n-l 



Here £'sa' ~ 10^^ ergs is the SN explosion energy and A'^ ~ 0.003 cm^-^ is the density of 
the hot phase of the ISM. A pulsar braking index of n = 2 ("LMC pulsar" type) will then 
increase the age of PSR B 1823- 13 to ~ 42 kyr, which would imply a size consistent 
with the predicted size of the unseen SNR shell of ~ 70 pc. 
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